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ABSTRACT 



An error probability analysis is done for a DFT based, M-ary frequency-shift 
keying (MFSK) communications system employing fast frequency- hopped spread 
spectrum signals. A linear combination procedure referred to as noise- normalization is 
employed at the receiver to minimize the effects of partial-band interference, which is 
modeled as additive Gaussian noise. The performance of the receiver is studied as a 
function of signal Doppler shift and type of windowing function used in the DFT. 

The use of fast frequency-hopped spread spectrum is found to improve the 
performance of the DFT based receiver in all but the most severe cases of Doppler shift. 
The use of a non-uniform window (i.e., a Hamming window) to improve receiver 
performance is effective only in the presence of large Doppler shifts. The amount of 
Doppler shift necessary to warrant the use of a non-uniform window depends on the 
amount of jamming noise power at the receiver, but is relatively insensitive to the 
frequency- hop rate used. In general, in the absence of any information concerning the 
nature of the received signal, a non-uniform window should be used because the 
performance degradation experienced at small Doppler shifts is insignificant compared 
to the performance enhancement gained at large Doppler shifts. 
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I . INTRODUCTION 

A. NON-UNIFORM WINDOWING OF THE DISCRETE FOURIER TRANSFORM 

The use of the discrete Fourier transform (DFT) to 
noncoherently detect signals employing M-ary frequency-shift 
keying (MFSK ) modulation has gained popularity recently due to 
advances made in high speed digital systems. At the current 
level of technology, however, the utility of the DFT based 
receiver remains limited due to the speed of real-time DFT 
algorithms. Even the performance enhancement obtained by using 
the fast Fourier transform (FFT) is not sufficient to overcome 
this limitation. However, certain applications - namely fast 
frequency-hopped spread spectrum systems - can effectively 
utilize the DFT receiver because lower data rates frequently 
correspond to higher immunity to narrowband interference. 
[Ref. 1 :p. 2014] 

Problems arise in using the DFT based receiver if the 
transmitted frequency does not match the assigned frequency of 
the DFT bin. This situation results when the signal frequency 
is different than that required by the detector, as when the 
signal experiences a Doppler shift prior to detection. This 
leads to a situation where transmitted signal energy falsely 
contributes to an incorrect DFT frequency sample: a phenomena 
known as self-induced crosstalk [Ref. l:p. 2015]. This problem 
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is not specific to the DFT based receiver. What is unique 
about this type of receiver is its ability to minimize self- 
induced crosstalk through the use of non-uniform windows. An 
example of a DFT based detector experiencing self-induced 
crosstalk is presented in Fig. 1.1. Figure 1.1(a) depicts the 
situation where the received frequency is perfectly matched to 
the first output bin of the detector, while Fig. 1.1(b) shows 
the results of slightly modifying the received signal's 
frequency. As can be seen, the magnitude of the bin 
corresponding to the signal is reduced, while the magnitudes 
of neighboring bins which represent other frequencies are 
significantly larger. Large amounts of frequency deviation can 
easily lead to unacceptable bit error rates. 

One method that is used to minimize the effects of self- 
induced crosstalk involves the use of nonuniform windowing; 
that is, time-sampled data is not equally weighted. The 
frequency response of typical nonuniform windows have broader 
main lobes than do rectangular windows (uniform weighting) but 
exhibit sharper rolloff characteristics. The sharp rolloff 
characteristics of nonuniform windows partially compensate for 
the effects of self-induced crosstalk, but in some cases the 
existence of a broader main lobe creates distortion where it 
would not otherwise exist [Ref. l:p. 2015]. Nonetheless, 
nonuniform windowing can be used successfully to combat large 
doppler shifts in low data rate environments as demonstrated 
in [ 1 ] and [ 2 ] . 
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(a) No Doppler (b) Doppler 



Figure 1.1: Outputs of a DFT (a) without a doppler shifted 
signal, and (b) with a doppler shifted signal. The adverse 
effects of a doppler shifted signal can severely degrade 
communications . 

B. FACTORS AFFECTING RECEIVER PERFORMANCE 

1. Fast Frequency-Hopping and Partial-Band Interference 

When the carrier of the transmitted frequency is 
changed in an apparently random manner during the transmission 
of a signal, the overall effect is to spread the energy of the 
signal over a large bandwidth. This spreading technique is 
called frequency-hopping and is illustrated in Fig. 1.2. Two 
types of frequency-hopping are considered: fast and slow 
frequency-hopping. Slow frequency-hopping systems change, or 
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hop, the carrier once per data symbol, while fast frequency- 
hopping systems hop the carrier more than once per data 
symbol. The higher the hopping rate, the more the frequency 
spectrum approximates that of true spread spectrum systems 
utilizing pseudonoise (PN) codes [Ref. 3: pp. 10-14]. 
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Figure 1.2: An example of hopping the carrier frequency of 
a signal. 



The intent of a hostile jammer is to disrupt 
communications. In order to effectively jam a frequency-hopped 
signal, a smart jammer may employ partial-band interference; 
although, this type of interference can also be caused by 
other unintended narrowband sources. In utilizing partial-band 
interference, the jammer concentrates all of the jamming 
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energy into a fraction of the total bandwidth and randomly 
hops this signal through the spectrum in the same fashion as 
the signal employing frequency-hopping. By doing this, the 
frequency hopped signal cannot always avoid the portion of the 
spectrum that is jammed; the jammer must be dealt with via the 
receiver . 

The assumptions involved in the analysis to follow are 
that the jammer is limited in power; that the interference can 
be modeled as additive Gaussian noise; and that the 
interference, when present, is present in all detection 
channels with probability y. Hence, the probability that 
interference is not present in any branch is 1-y. Thus, y 
represents the fraction of the spread bandwidth that is 
jammed. If the partial-band interference is assumed to have a 
power spectral density of N^/2 over the entire spread 
bandwidth, the amount of interference present at the receiver 
due to partial-band interference is y -1 Nj/2 with probability y 
and 0 with probability 1-y. [Ref. 4:p. 3] 

In addition to partial-band interference, the system is 
assumed to be corrupted by additive white Gaussian noise with 
a power spectral density of N Q /2. Accordingly, the power 
spectral density of the total noise at the receiver input is 
N 0 /2+y _1 Nj/2 with probability y and N Q /2 with probability 1-y. 
If the receiver equivalent noise bandwidth is W Hz, then the 
total noise power at the receiver is N Q W with probability 1-y 
and is (y -1 Nj+N 0 )W with probability y. [Ref. 4:p. 3] 
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2 . Fading 

In addition to partial-band interference, the 
signal is assumed to be further degraded by fading. Fading is 
caused by components of the signal arriving over different 
paths which combine to form the total detected signal. Due to 
the different path lengths traversed by the various signal 
components, the combined signal amplitude may be greater than 
or less than the expected amplitude without fading due to 
constructive and destructive interference. Fading severely 
degrades the performance of a real world communications system 
and must be taken into account in any plausible receiver 
design and analysis. 

For this analysis, the signal is assumed to experience 
slow, frequency non-selective, Rician fading. Slow fading 
refers to the fact that the amplitude of the signal is assumed 
to remain constant at least over the duration of one hop, 
while frequency non-selective implies that all frequency 
components of the signal experience identical fading. Rician 
fading is characterized as having a diffuse signal component 
and a direct signal component. Assumme that the amplitude of 
the received signal in a Rician fading channel is a/2 . Then 
the probabilty density function of the random variable A is 
given by 



fAa)=- 



2of 



i A 



aT 



) u(a) 



( 1 ) 
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where u(-) is the unit step function, r 2 is the power in the 
direct component of the signal, 2a f 2 is the power in the 
diffuse component, and I Q ( • ) is the modified Bessel function 
of order zero. As can be seen, (1) reduces to the well known 
Rayleigh distribution when the value of the direct component 
is zero. [Ref. 5] 

C. THE NOISE-NORMALIZED DETECTOR 

A noise-normalized detector is a linear receiver that uses 
the received noise power to normalize the output of each 
branch before a decision is made regarding which signal is 
present. This type of receiver is combined with fast frequency 
hopping to improve overall system immunity to partial-band 
interference. The noise-normalization is implemented before 
the hop receptions are combined to form the detection 
statistics. This type of detection scheme tends to lessen the 
influence on the overall decision statistic of each hop when 
interference is present, while increasing the influence of 
each hop when no interference exists [Ref. 4:p. 2]. Because 
the output statistics of each branch of the receiver are 
independent of one another, closed form solutions for the 
receiver performance can be obtained. 

This thesis studies the effects of nonuniform 
windowing on a noise-normalized DFT based receiver where the 
signal experiences partial-band interference, Doppler shift, 
and Rician fading. The effects of fast-frequency hopping the 
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carrier frequency are analyzed and the results for various hop 
rates discussed. As in [1] and [2], a statistical analysis is 
done to determine the amount of frequency deviation necessary 
to warrant the use of a nonuniform window over a rectangular 
window. 
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II. ANALYSIS OF THE NOISE-NORMALIZED DFT BASED DETECTOR 



The receiver under consideration is shown in Fig. 2.1. The 
incoming signal hop is first brought back to baseband and then 
sampled at the Nyquist rate. The DFT of the signal is computed 
using the appropriate windowing function and the magnitude of 
each output bin is squared. Hence, in the absence of both a 
Doppler shift and a nonuniform window, this receiver is 
functionally equivalent to a conventional MFSK receiver 
employing quadratic detectors. At the same time, the total 
noise power at the receiver is measured, modified 
appropriately by the windowing function of choice, and 
inverted to form the noise-normalization term. Each bin output 
is multiplied by the noise-normalization term. The decision 
statistics are then obtained by summing over all hops 
comprising the transmitted symbol. Finally, the detection 
statistics are compared and the largest one used to determine 
the estimate of the symbol sent. Errors result when the 
detection statistic in a channel with no signal is greater 
than the detection statistic of the signal channel. A detailed 
analysis of each part of this design follows. 
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Figure 2.1: The DFT based noise-normalized receiver 
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A. DETERMINATION OF THE PROBABILITY DENSITY FUNCTION FOR EACH 



DFT OUTPUT BIN 

The signal input to the receiver in Fig. 2.1 every T h 
seconds is of the form 



where CKt_<T h , a is a Rician random variable having the 
probability density function of (1), aV2 is the Rician 
distributed signal amplitude, f c is the carrier frequency 
hopped every T h seconds, f m is the intended symbol frequency, 
f is a Doppler shift, and n'(t) is zero mean, additive 
Gaussian noise. The signal frequencies are chosen to be 
orthogonal during the time interval 0<t<T h and are obtained 
using the relation 



where m'=l,2,..,M and is an integer representing the 
spacing between each of the M signals. To simplify the signal 
representation, the following substitutions are made: 



r ( t) =aj 2cos [2 k ( f c +f m +f ') t] +n'( t) , 



( 2 ) 



1 A £ (m'-l) 

r m /= + 



(3) 




g£l+A f (m'-l) 



(4) 



and 




(5) 



so that after the low-pass filter (2) becomes 
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r'( t) =a/ 2cos [ 



] +J2 ( t) , 



( 6 ) 



2n (q+e) t 
T h 

where e is the fraction of Doppler shift with respect to the 
bandwidth of the baseband information signal, q is an integer 
related to the signal frequency f m . , and n(t) is the filtered 
noise. The discrete form of the received signal is obtained by 
using a sampling rate of N/T h , where N is the number of 
samples taken during the symbol reception: 

rip) =aj2cos[ 27t( y e) ^ ] + n(p) , (7) 

N 

where p=0,l,..,N-l and is the discrete signal sample. [Ref. 
2 : pp . 7 - 9 ] 

The output of each bin is a complex quantity and can be 
separated into real and imaginary parts: 

Yjk\a) =Y Rm (k\a) +jY Im (k\a) , ( 8 ) 

where k, which ranges from 0 to N- 1, is an integer 
representing the DFT output bin and F m (k|a), y Rm (k|a), and 
7 Im (k|a) are respectively the DFT output, the real part of the 
DFT output, and the imaginary part of the DFT output for 
channel k and frequency hop m conditioned on the Rician 
distributed random variable a. The real and imaginary parts of 
the DFT are calculated separately 



12 



( 9 ) 



=Y^ 0 r{p) wip)c os 

=V*™ [v^acos[ 2 *(g +e >P ] +/3(p)] w(p)cos[-^ffi] , 

■‘-'P* o AT AT 

and 



y im(*|a) =££J * <P> *(P> sin[-^£] (10) 

-JTi t/2acos[ 27t( ? - t - e . ) - P ] +n(p)] w(p)sin[^£] , 
i—tp-a at n 

where w(p) represents the windowing function. Assuming slow 
Rician fading so that the amplitude of the signal remains 
constant during each hop interval, the conditional probability 
density function for the random variable X km is 



-<**1 



A^ X km 






-mL.) 



2o, 






( 11 ) 



where X km is the magnitude of bin k and hop m given by 



*ta=|r a (ic| a) | 



( 12 ) 



p km is the magnitude of the expected value of bin k and 
freguency hop m given by 



M'fan l^cml 

and a x 2 is the variance of bin k given by 

<* x 2 = \z[\Yjk\a) -Y m {k\a) \ 2 ] . 

[Ref. 1 :p. 2016] 



(13) 



(14) 
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The mean of X km for a fixed amplitude a is easily 
determined by observing that the only random quantity is n(p) 
which is zero-mean. This leads to the conclusion that the mean 
is simply the output of the DFT without noise. The real and 
imaginary parts of the mean conditioned on a are 



L Rm 



net aT-a#££, ^ c ° s < ) cos ( ) 



( 15 ) 



tun- 



Rm 



and 






N 



N 



( 16 ) 






Im 



An expression for p km 2 conditioned on a is then determined by 
substituting (15) and (16) into (13) and squaring both sides: 



Vkm 2= ™R + ™l 

=a 2 (m' 2 R +m'\) (I 7 ) 

where /3 km is the magnitude squared of the mean of the DFT 
output without fading when the signal power is normalized to 
unity, that is, 



e *.' 2 < ££><*>■ cos < **Bte*L ) cos ( - 22 fi* ) J 3 
.[V *- 1 cos( 2 l, P<q^ )sin(^£jl)] 2 ) . 

N N 



( 18 ) 



[Ref. 2:pp. 10-11] 

The variance of X km is obtained by substituting (9), (10), 
(15), and (16) into (14) and simplifying 
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( 19 ) 



^ x 2 =\ E ^\ Y ^k\a) -T^JTc\a) +j (Y Iln (k\a) -Y Im (k\a) ) | 2 ] 

=E[ lE^o "(P)n(p)c os(-^£) 

+j Y^> wip) n {p) sin ( ^jr~ ] l 2] • 

The magnitude squared of (19) becomes 

°x = \ E[ "(P)n(p) cosi^S .)) 2 

+ (IT' 1 w(p) n ( p ) sin ( -^£2 ) ) 2 ] 

*-^p - o TV 

=B[ 52*-; V*'! »'<P> "'<«> > cos < ■ 22 ^ 2 )n(p)n(s) 

x-^p= o tss=o n N 

+IT 'l IT 'l w(p) w(s) sin ( 2Tt f- P ) sin( 2u ^ s ) n(p) n(s) ] , 
z^p=o o N N 

where s is an integer having the same range of values as p. 
Because the expected value operator is a linear operator and 
the only random quantities are the noise functions, (20) 
simplifies to 



o x 2 =ir 1 y N 1 w(p) w(s) cos ( ^luJsp ) 

x 2 Y~>p=o Z^/s-o ^ i\T 

• cos (^EM) E[n(p) n(s) ] 

+ — y?"" 1 IT' 1 Wp) v(s) sin <_2«icp ) 

2 ^p=o ^s = o ^ N 

. sin (-^£)E[n(p)n(s)], 

N 



( 21 ) 



The autocorrelation function is defined as 



RrmiP.s) =E[n(p)n(s)] , ( 22 ) 

which, because the noise functions represent white noise, is 

R nn {p, s) =WN c b (p-s) , ( 23 ) 

where 6(-) is the Dirac delta function, W represents the 
equivalent noise bandwidth of the low-pass filter at the front 
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of the receiver, and N t is the power spectral density of the 
total noise seen at the receiver [Ref. 4]. The quantity N t 
equals either N Q when jamming is not present or N Q +N^/y when 
jamming is present. The expected value needed to solve (21) is 
obtained by substituting (23) into (22): 



The final expression for the variance ' of X km is obtained 
by substituting (24) into (21) and simplifying: 



The variance in each bin is affected by both the power 
spectral density of the external noise and the choice of 
windowing function used in the DFT. 

The probability density function for the output of each 
bin is found by substituting (17) into (11) and integrating 
over the range of the Rician distributed amplitude: 



E[n (p) n is) ] =WN t 6 (p-s) . 



( 24 ) 
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f xJ x kn) =f~ f x km \A ixJ\a)f A {a)da 



(x^a 2 p to ) 



-r-s i 

Jo o / 



2 <Jv 



- (a 2 +I' 2 ) 



V' 



ax ; 



JcnVP jcm 



- 5 -e 2 «' ! J 0 (^)da 



( 26 ) 



v icm 



O^O 



2 rr 2 






.2(J5= + J1) 

2 °* 2 O f 2 



•fae ' 2o * 2 20f2 ’ J 0 ( ax ta ^ ) J 0 ( ) da . 

Jo (i 2 n - 



By making the substitutions I Q (x)=J Q ( jx) , and by defining 



q2 a fikn + i_ 



20/ 20/ 

_ P*m°f 2+g k 2 
2o/o/ 



(27) 



and using the integral (from [7]) 



f xe~° 2t2 J^ {ax) J v (bx) dx=— ^—e 4£?2 , 

Jo 2 Q 2 2 Q 2 

the probability density function is 



(28) 






x 



km 



xZftf 2 '**** 2 



2o 



2Q 2 o x 2 o f 2 






( — ) ( 

iQ 2 



Xkjr$k/n 



o f 4 j / -^ion^V P tan 

0 2 Q 2 o 2 o f 2 



) • 



(29) 



The next step involves substituting (27) into (29) which gives 
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( 30 ) 



f xjx jJ=( 



X 



-^LV-fox 2 ^ xLv 



km 



) e 



2a x°t 2a x‘ °f Wkm°t : +a x 1 



Pjcm 0 / +0 x 2 

, 7 - / ^knP yPjcm \ 

V ft 2 2 * 

P*m°f + °* 



Equation (30) simplifies to 



1 / an * P Arirf - , 

p 2 Vkm°f 2 +°x 2 T ( X knF'JPkm \ 

.X*7n v ^km* q o 2 "^^'ft'2 2^ 

PtaOj *°* P»,°/*°x 






X; 



fan 



( 31 ) 



In order to clarify the notation, the following substitutions 
are made: 



a km 2 *Vkm a f 2+a x 2 ' 



( 32 ) 



and 

“*.* 1715 ^. ( 33 ) 

With these substitutions, (31) is seen to be Rician 
distributed 



f xJ x J»>- 






km 



2q v 



9 km 



V 



Xkm a km 



’km 



) , 



( 34 ) 



and reduces to the result obtained in [1] for the case of a 
nonfading channel (2o f 2 =0 and r=l). [Ref. 2:pp. 11-12] 
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B. ANALYSIS OF THE NOISE-NORMALIZED RECEIVER 



1. Determination of the Probability Density Function for 
the Detection Statistic 

As illustrated in Fig. 2.1, the random variable Z km is 

given by 



Z km 




( 35 ) 



where £ km is the random variable for bin k and frequency hop 
m. Because all the random variables X km are independent, the 
transformation of random variables is accomplished by using 
the relation (from [8]) 



■ 2*1 ( Z *J 



0 ) ) 

Z'^kmiO) ) 



( 36 ) 



where 



Z km=£r( x km> > ( 3? ) 

g'(x km ) is the derivative of g(x km ) with respext to x km , and 
x km(0) t ^ le rea l root of (35). The denominator in (36) is 
found by substituting (35) into (37) and differentiating: 
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a' lx ) - 2X km(°) 
3 ' x km ( 0 ) ' 





o 2 

X 


2 


^ km 




°x 2 



( 38 ) 



The transformation of random variables for the noise- 
normalized detection statistic is determined by substituting 
( 38 ) into ( 36 ) : 



^Zkm ( z kn ) 2 \ 






( 39 ) 



■‘km 



The probability density function of Z km is obtained by 
substituting (34) into (39): 



- l°x 2z km+* 2 kJ 



^ 7. ( ^Irrr) 



2 



2°L 



*n< 






( 40 ) 



, km 



The actual detection statistic needed for the 
comparison is Z k , which is defined as 

Z k = Y! m = 1 Z km’ < 41 ) 

where L is the number of frequency hops per symbol. Since some 
hops are jammed while others are not, the probability density 
function of Z k conditioned on the number of jammed hops is 
simply the multiple convolution of the probability density 
functions of the jammed and unjammed hops. For instance, if a 
symbol has three hops and two of them are jammed, the 
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probability density function of the detection statistic is 
determined as follows: 




where <g> represents convolution. In general, the probability 
density function of Z k conditioned on 1 out of L hops jammed 
is 



where ®c n represents a c n -fold convolution. [Ref. 4:p. 7] 

2. Determining the Probability of Bit Error 



accomplished, the next step involves calculating the 
probability of bit error. An error will occur if the power in 
a channel other than the signal channel is greater than the 
power in the signal channel. Mathematically, this is expressed 
as 



where i,j=l,...,M and i*j ; and the signal is assumed to be in 
branch i.[Ref. 4: p. 7]. The double integral in (44) 
determines the probability that the power present in the 
signal channel is less than the power present in a noise 
channel . 




Once the random variable transformation is 




(44) 
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Because the distribution of power in each noise 
channel is not identical when Doppler shift is present, the 
comparison in (44) must be carried out for each noise channel 
and the results combined to give the total probability of 
error. An upper bound to the total conditional probability of 
error is found using the union bound: 

Ps u^Yj 3 = \ PriZ^Zj) . ■ ( 45 ) 

j*i 

where Ps iL is the probability of symbol error given 1 hops are 
jammed with the signal present in channel i. [Ref. l:p. 2017] 
Because the distribution of power in each bin varies 
with the location of the signal, the amount of Doppler, and 
the relative direction of the freguency shift, an average must 
be taken over all signal channels and direction of Doppler 
shift : 



Ps, (e) x Ps n (+e) +Ps n <" e > ' < 46 > 

where Ps L (e) is the total probability of symbol error given 1 
hops are jammed conditioned on the signal experiencing a 
Doppler shift of e and Ps u (*) is the probability of symbol 
error given 1 hops are jammed with the signal present in 
channel i conditioned on the signal experiencing a Doppler 
shift of (+£) or (-£). The union bound on the conditional 
probability of symbol error is calculated by combining (45) 
and ( 46 ) : 
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( 47 ) 




The total conditional probability of symbol error assuming 1 
hops are jammed is calculated by combining (44) with (47): 



[Ref. l:p. 2017] 

The total probability of symbol error for the noise- 
normalizing receiver in the presence of partial-band 
interference is 



where Ps(e) is the total probability of symbol error [Ref. 
4:p. 6]. To determine the worst case probability of symbol 
error, (49) is evaluated for different values of y, the 
results compared, and the highest conditional probability of 
symbol error recorded. The worst case probability of bit error 
is : 



where Pb(e) is the probability of bit error. 

An exact solution for the conditional probability of 
symbol error is calculated by evaluating the probability of 
not making an error (i.e., the probability that the power in 





( 49 ) 




( 50 ) 
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all noise channels is less than the power in a signal 
channel ) : 




( 51 ) 



The total conditional probability of not making an error is 
found by the following relation: 



where Ps i [ ( 1 ) is the probability of not making an error given 
that the signal in channel i with 1 hops jammed and the signal 
experiencing a Doppler shift of e. The probability of making 
an error given that the signal is in channel i with 1 hops 
jammed is then given by 



which is then used in (46), (47), (48), (49), and (50) to 
obtain the exact probability of bit error. Due to the enormous 
amount of computing power needed to numerically integrate (51) 
and the desire to compare the results with those obtained in 
[1] and [2], the conditional probability of error based on the 
union bound was chosen for this thesis. 



Ps i2 <e> =n>i Pr(Zj<Zi\e) 



( 52 ) 




Ps n (e) =l-Ps i2 (e) . 



( 53 ) 
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III. NUMERICAL PROCEDURE 



A. PARAMETERIZATION OF VARIABLES 

In order to compare the performance of the noise- 
normalizing receiver with a conventional DFT/MFSK receiver, 
the number of points used in the DFT, the cut-off frequency of 
the low pass filter, and the spacing between the M signals are 
chosen to be consistent with those used in [1] and [2], The 
following definitions result from adopting these conventions: 



N&4M, 


(54) 


W*2M, 


(55) 



and 



A f A2, (56) 

where N is the number of signal samples per hop, M is the 
modulation order, VJ is the effective noise bandwidth, and 
is an integer representing the message frequency spacing. 

A hop time, T h , of one second results from sampling the 
received signal at the Nyquist rate. When this value of T h and 
(56) are substituted into (3), it is determined that f m 
consists only of odd values ranging from 1 to 2M-1. 
Furthermore, the spacing of frequencies in the DFT is 
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A f= 



( 57 ) 



£s 

N 
AM 
AM 

= 1HZ, 

where Af is the frequency resolution of the DFT [Ref. 9]. 
Because of these choices one empty bin exists between each 
possible DFT signal output. 

A parameter of concern is the energy-per-bit to thermal 
noise power spectral density ratio defined as E b /N Q , where E b 
is the transmitted energy per bit. The signal-to-noise ratio, 
which is a measure of the ratio of signal power to noise 
power, is defined as 

r 2 +2o f 2 /cox 

SNR& (58) 

N 0 W 

For an M-ary system, symbols are sent instead of bits. The 
power per symbol is the sum of the power in the direct 
component and the power in the diffuse component: 

E s = (2a f 2 +r 2 ) T s , (59) 

where E s is the transmitted energy per symbol and T $ is the 
symbol duration time related to the hop duration time by 

T s =LT h . (60) 

The energy-per-symbol is related to the energy-per-bit by 
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E s =E b log 2 (M) . ( 61 ) 

All the variables used in this thesis are expressed as per-bit 
quantities; therefore, a relationship between the energy-per- 
hop and the energy-per-bit is needed. Because each symbol is 
comprised of L bits, 



E h = 




( 62 ) 



where E b is the energy-per-hop. A relationship between the 
energy-per-bit to noise power spectral density ratio and the 
signal-to-noise ratio is obtained by combining (58) with (59) 



£b 

^0 



WT S 

log 2 (M) 



SNR. 



( 63 ) 



Proper choice of the symbol duration time (i.e., T s =l) and the 
equivalent noise bandwidth (i.e., W=log 2 (M)) gives E b /N Q =SNR. 
In general, however, these two ratios are not equal and must 
be kept separate in the receiver analysis. 

An expression relating the energy-per-bit to noise power 
spectral density ratio to the signal-to-noise ratio is 
obtained by combining (58), (60), and (63) and simplifying: 

E b _ M (L) (T h ) r 2 +2o f 2 (64) 

N 0 log 2 (M) N 0 W 

For the case of a receiver operating in a jammed environment, 
N q is replaced with N Q +N^/y when the hop is jammed and (64) 
becomes 
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( 65 ) 




r 2 +2o/ <R0 (D (T h ) 

W(N 0 +Nj/ y) log 2 (M) 



The relationship between the energy-per-bit to noise power 
spectral density and the signal-to-noise ratio is a function 
of the effective noise bandwidth of the system, the number of 
freguency hops per symbol, the hop duration time, and the 
level of M-ary signalling used. 

Another parameter for analysis is the direct to diffuse 
ratio 



The following expressions for the direct and the diffuse 
signal power components result from combining (65) and (66) 
and simplifying: 



In doing this analysis, the values of E b /N Q , N^/N QI y, W, T h , 
and DTD are assumed and used to determine the values of T 2 and 



DTD& 




( 66 ) 



f 2(l+DTD)T h 



( 67 ) 



and 



p _ ffh / DTD_ \ 
’ T h 1 +DTD 



( 68 ) 
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B. CONVOLUTING THE PROBABILITY DENSITY FUNCTIONS 

In the case of fast frequency-hopping, the probability 
density function of the total detection statistic is the L- 
fold convolution of the probability density functions of the 
detection statistic of each hop. As shown in (49), the total 
probability of error represents all possible combinations of 
jammed and unjammed hops. Based on this, a total of 2M(L+1) L- 
fold convolutions are required for each possible signal branch 
and for each combination of jamming power and fraction of 
bandwidth jammed that are studied. These convolutions are not 
required in the case of slow hopping because the total 
detection statistic is equal to the detection statistic of the 
single hop. These multiple convolutions are generated by 
taking advantage of the properties of the Laplace transform. 

The Laplace transform of the probability density function 
of Z im is 

i m 



where n=l,2 and n = 1 represents the condition of a jammed hop 
and n=2 represents the condition of an un jammed hop. The 
Laplace transform of interest is obtained by substituting (40) 
into (69) 



(Z^)] =F z (n) (S) 



( 69 ) 
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( 70 ) 



,<n> = f 
'in J 0 



- ( 0 x 2z in > +0 jjn 2<nl > 



o 2a- 2(n> 
^ iw 



20 ™ <n> J n ( a ^ < ” >q ^V " Ja ) e " SZira dzi 






Tn) 



szjrt (fl) 



OilB 



2 (n) 



2 (n) 



2 o J „ 2 <"» r- - 



2o 



2 (n) 



/• 









,2 ( 73 ) 



" )dzi?. 



Equation (70) is evaluated by defining 



£= 



ja, J n) o, 



* j/n x 



® im 



2 (n) 



(71) 



and 



C 2 =s+ 



2o 



2(n) ' 



(72) 



and making the substitutions u 2 =z im and I 0 (x) =J 0 ( jx) , and using 
the integral (from [10]) 



/ 

Jo 



e~ c2tj2 uJ 0 (Bu) du = — — e 4 ° 2 
0 2C 2 



The Laplace transform of interest is 



F? (n) ( s ) = 



A< n> 



A (n> 



where 



(73) 



(74) 



A (n> =2o ijn 2 !n> S+o x 2 (75) 

[Ref. 4:p. 6], The probability density function of Z- is 

determined by using the fact that convolution in the time 
domain translates to multiplication in the s domain: 
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( 76 ) 



The probability density function of is given by 

f z Az i )=g-H(F z m(s)) 1 (F zl 2j(s)) L - 1 ] , ( 77 ) 

where L is the number of hops per symbol, 1 is the number of 
jammed hops, and Sf 1 is the inverse Laplace transform 
operation. In general, (77) cannot be evaluated in closed form 
and numerical subroutines are used to determine the inverse 
Laplace transform at each point. 

For the case of all hops jammed or unjammed, (77) reduces 
to 



£ z .„(z 1 (j,, )=se- 1 [(^. to ,(s)) i ] 

= [~(F z <n> (s)) L e sz <"ds. 

Jo 



( 78 ) 



The function in the s domain in (78) is obtained by raising 
(74) to the L power 



, (n) 



Ls 



[ F 7 in) (S) ] L =( 



2 a ‘ 

im 



I (n) 



) L e 



2 <J im 2 



( D ) 



s+o. 



Equation (79) is factored as 



( 79 ) 
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( 80 ) 



La,- 



2 (n) 



[F„<„>(s)] L =( 



O _2 (n) 

2o im 



v r 2o • 2 (n) 

) L e 



r ..2{n) /« 2 / / „ „4(n) s 
La ii» a x/( 4a ii7» ) 



(s+- 



2o i« 



2(n) 



)‘ L € 



2^ 2(D> 



The needed inverse Laplace transform is obtained by combining 
the transform property that 

if* 1 [F(s+a) ] =e' az f (z) , ( 81 ) 

the inverse Laplace transform relation (from [11]) 

S£ -1 [ ( — ) L e =(-^) "^ i J i . 1 (2 v /Si) , < 82 ) 

s b 



with the substitutions 



a=- 



2o 



2(n) 



and 



b= 



La . 



2 (n), 






4 (n) 




(83) 



(84) 



After simplifying, the resulting closed form expression for 
this special case is 



f z w (zj n) )=- 



fl (n) 

2o. 



r-i 

) 2 

2 (n) 



2o Jm 2(D) 



(z+Q In> ) 






2 (n) 



VzQ 17 ^) , 



(85) 



where I L _ 1 is the modified Bessel function of order L-l and 
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( 86 ) 



„2(n> 

® int 



[Ref. 4 : p . 7]. 

C. NUMERICAL INTEGRATION TECHNIQUES 

A double numerical integration is needed to evaluate (48), 
which can be very time consuming. By taking advantage of the 
properties of the Laplace transform, the double numerical 
integration is reduced to a single numerical integration. The 
amount of computing time needed is greatly reduced since 



and (76) can be divided by s as needed to carry out the inner 
numerical integration. 

The use of the integration technique illustrated in (87) 
is the primary reason a union bound solution is calculated 
instead of an exact solution. As seen in (44), the union bound 
solution requires implementation of (87) only once. The exact 
solution, as seen in (51), requires using (87) M-l times 
since, for the DFT based receiver, the bins that ideally 
contain no signal will in general have different probability 
density functions due to Doppler shift and non-uniform 
windowing. Although using (87) to integrate the probability 
density function significantly reduces computing time compared 
to direct numerical integration, the amount of computing 




( 87 ) 
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resources needed to calculate the exact solution is much more 
than that needed to determine the union bound solution. The 
fact that only a limited amount of these resources was 
available to generate results justified the decision to 
evaluate the union bound solution. 

The remaining integration is accomplished using Simpson's 
rule, which is chosen because of its accuracy as compared with 
the trapezoidal rule. Figure 3.1 shows how Simpson's rule is 
applied to integration. By dividing the z-axis into equally 
spaced elements and evaluating the function at each of these 
points a close approximation of the integral can be obtained 
through the relation 

j *' "f(z ) dz=A^- (f(z 1 )+4f(z 2 )+2f{z 3 )+4f(z 4 )+2f(z 5 ) 

+ . . . +4f{z N _ 1 ) +f[z N ) ) 

where an odd number of points must be used [Ref. 12 :p. 350]. 
By choosing an endpoint, varying the spacing between points, 
and repeatedly evaluating the integral until no difference is 
observed between successive iterations, an optimal Az is 
obtained. A close approximation to the semi-infinite integral 
is obtained by using that value of A z while increasing the 
value of the upper limit on the integral until no difference 
is detected between successive iterations. 
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Figure 3.1: Example of dividing a function into equally 
spaced segments suitable for use with Simpson's rule. 
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IV. RESULTS 



A. PREVIOUS RESULTS 

Studies concerning the effects of non-uniform windowing on 
a DFT based receiver in the presence of a Doppler shifted 
signal have been performed and documented in [1] and [2]. Each 
study draws conclusions concerning the effectiveness of the 
DFT based receiver under different circumstances. The results 
of these studies are included in order to better understand 
the conclusions of this thesis. 

Reference [1] studies the effects of non-uniform windowing 
on a DFT based receiver in a nonfading channel with a Doppler 
shifted signal, where the probability of bit error as a 
function of the signal-to-noise ratio is obtained. It is 
demonstrated that using a non-uniform window (specifically, a 
Hamming window) results in a degradation of performance as 
compared with the performance of a uniform window for the 
cases of little of no Doppler shift on the incoming signal. 
However, performance enhancement is obtained as Doppler shift 
increases. For the case of M=8, a Doppler shift of e=.33 is 
sufficient to warrant the use of non-uniform windowing. 
Furthermore, the performance improvement increases as the 
amount of Doppler shift in the signal increases. [Ref. 1: pp. 
2017-2018] 
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The effects of a Rician fading channel on the DFT based 
receiver are investigated in [2], Another parameter for 
analysis is the direct-to-dif f use ratio discussed in the 
previous chapter. In a manner similar to that developed in 
[1], data is generated relating the probability of bit error 
to the signal-to-noise ratio for different direct-to-dif fuse 
ratios. For large direct-to-dif fuse ratios (i.e., conditions 
corresponding to those studied in [1]), the results correlate 
well. However, as the direct-to-dif fuse ratio is decreased, 
the non-uniform window gives better performance than a 
rectangular window only when the Doppler shift decreases. The 
presence of a strong diffuse signal, as illustrated by a small 
direct-to-dif fuse ratio, limits the benefits of using a 
rectangular window. 

B. EFFECTS OF FREQUENCY-HOPPING AND PARTIAL-BAND INTERFERENCE 

Graphs of the worst case probability of bit error as a 
function of signal-to- jamming noise ratio for fixed signal-to- 
therraal noise ratios are constructed for one, two, three, and 
four hops per symbol. These results are obtained by 
numerically evaluating (40), (48), and (49) for values of y 
ranging from 1.0 to 0.001 and retaining the worst case 
performance for each value of F b /N- . These results are plotted 
in Figs. 4.1 through 4.4. The signal-to-noise ratio, direct to 
diffuse ratio, and the modulation order are chosen to enable 
direct comparisons with the results obtained in [4] which 
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Figure 4.1: Worst case probability of bit error versus 
signal-to- jamming noise ratio for L=1 , M=4, direct-to- 
diffuse ratio=10, and SNR=13.35 dB . 



38 



Eb/Nj (dB) 




Figure 4.2: Worst case probability of bit error versus 
signal-to- jamming noise ratio for L=2, M=4, direct-to- 
diffuse ratio=10, and SNR=13.35 dB. 
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Figure 4.3: Worst case probability of bit error versus 
signal-to- jamming noise ratio for L=3, M=4, direct-to- 
diffuse ratio=10, and SNR=13.35 dB. 
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Figure 4.4: Worst case probability of bit error versus 
signal-to- jamming noise ratio for L=4, M=4, direct-to- 
diffuse ratio=10, and SNR=13.35 dB . 
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analyzes a conventional (non-DFT) noise-normalized MFSK 
receiver with quadratic detection. The results obtained here 
for the case of no Doppler shift (£=0) and a rectangular 
window should be similar to those obtained in [4]. 

Differences in the results obtained in [4] and those 
obtained in this study are primarily due to the facts that the 
equivalent system bandwidths used in [4] and this thesis are 
different and the results in [4] are exact while the results 
in this study represent a union bound. The equivalent system 
bandwidth used in [4] is chosen to be minimal while the 
equivalent system bandwidth used in this work is chosen to 
enable increased frequency spacing between the DFT output bins 
as well as to allow comparisons to be made with the results in 
[1] and [2]. From this, it is expected that the data generated 
for this thesis will predict slightly worse performance than 
that predicted in [4]. This is observed in all cases. 

Just as reported in [4], the performance of the DFT based 
communications system is enhanced as the number of frequency- 
hops per symbol is increased. This is especially true for the 
cases of little or no Doppler shift. For cases of large 
Doppler shift (i.e., £=0.5), the performance of the receiver 
is slightly poorer as the hop rate is increased when non- 
uniform windowing is used. The error introduced by self- 
induced crosstalk overcomes the benefits usually associated 
with fast frequency-hopping. This is clearly seen in Fig. 4.1 
and Fig. 4.4 for the case of a rectangular window and e=0.5. 
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For E b /N- = 41.0 dB, Pb= 4.41xl0 -2 for a hop rate of one hop per 
symbol, and Pb=6.65xl0' 2 for a hop rate of four hops per 
symbol. The use of a Hamming window in this situation improves 
receiver performance as the number of hops per symbol 
increases, but the benefits are small compared to the 
complexities needed to implement a fast frequency-hopped 
receiver . 

For the case of slow frequency-hopping, the results in 
Fig. 4.1 approach those reported in [2] as E b /Nj becomes very 
large. This is expected because the decision statistic 
consists only of one component due to the fact that the 
carrier frequency is changed once per symbol. Furthermore, 
receiver performance degrades as either the jamming power or 
the Doppler shift increases. As reported in [1] and [2], the 
performance of the Hamming window is worse than that for a 
rectangular window for small Doppler shifts, but outperforms 
the rectangular window for large Doppler shifts. The value of 
e at which the use of a Hamming window gives better 
performance than the use of a rectangular window (the 
crossover fractional frequency offset) is a function of the 
amount of jamming power at the receiver. At E b /W.=1.0 dB, for 
example, with a rectangular window, Pb=0.413 with £=0 and 
Pb=0.687 with £=0.5; while, with a Hamming window, Pb=0.524 
with £=0 and Pb=0.668 with e=0.5. However, at E b /N ^=4 1.0 dB, 
with a rectangular window, Pb=4.59xl0 -4 with £=0 and 
Pb=4.41xl0' 2 with £=0.5; while, with a Hamming window, 
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Pb= 1.26xl0 -3 with e=0 and Pb=7.87xl0 _3 with e=0.5. For large 
Doppler shifts, the use of a Hamming window clearly improves 
receiver performance over a broad range of E b /N j . 

This same phenomena exists for the cases of two, three, 
and four hops per symbol. In these cases, however, the larger 
hop rates have little effect on the value of the crossover 
fractional frequency offset. This is illustrated in Figs. 4.5 
through 4.8 which plot the worst case probability of bit error 
as a function of the fractional frequency offset for 
E b /Nj= 16.0 dB and L=l,2,3 and 4. These results are summarized 
in Fig. 4.9, which plots the crossover fractional frequency 
offset as a function of the hop rate for various values of 
£ b /N- . As can be seen, the values of the crossover fractional 
frequency offset remain in a small band for each value of 
E b /Nj , and increase as the amount of jamming power at the 
receiver increases. 

From Fig. 4.9, it can be seen that the crossover 
fractional frequency offset does not follow a smooth curve as 
the hop rate is varied. In viewing these curves, it must be 
remembered that the results in Figs. 4.1 through 4.4 represent 
the worst case performances of union bound solutions. These 
worst case results are obtained by varying the value of the 
fraction of spectrum jammed (y) and retaining only the worst 
case performance. The value of y giving worst case performance 
at one hop rate is not generally the same value at a different 
hop rate. The fact that only a representative number of 
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discrete values of y are used to approximate the continuous 
range of possible y's contributes to this situation. 
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Figure 4.5: Worst case probability of bit error versus 
fractional frequency offset for the case of L=l, M=4, 
SNR=13 . 35 dB, and E b /N j =16.0 dB. 
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Figure 4.6: Worst case probability of bit error versus 

fractional frequency offset for the case of L=2, M=4, 

SNR=13 . 35 dB, and E b /N j =16.0 dB . 
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Figure 4.7: Worst case probability of bit error versus 
fractional frequency offset for the case of L=3, M=4, 

SNR=1 3 . 35 dB, and E b /N j =16.0 dB. 
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Figure 4.8: Worst case probability of bit error versus 

fractional frequency offset for the case of L=4, M=4, 

SNR=1 3 . 35 dB , and E b /Nj=16.0 dB . 
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Figure 4.9: Crossover fractional frequency offset versus 
frequency-hops per symbol for the case of M=4, and 
SNR=13 . 35 dB. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



As the frequency-hop rate of the DFT based communications 
system with noise-normalization increases, receiver 
performance improves in the absence of large Doppler shifts. 
However, when large Doppler shifts exist, the errors 
introduced as a result of self-induced crosstalk counteract 
any benefits gained from the use of fast frequency-hopping. 
The advantages gained by using a non-uniform window in the 
presence of large Doppler shifts improve the situation, but 
not enough to warrant the added hardware complexities inherent 
to the application of fast frequency-hopped spread spectrum. 

The amount of Doppler shift necessary to warrant the use 
of a non-uniform window is primarily a function of the jamming 
noise power spectral density. Holding all else constant, as N ^ 
increases, the amount of Doppler shift necessary to justify 
the use of a non-uniform window increases. At very small 
values of E b /N ■ the use of a non-uniform window results in 
degraded performance for all but the most extreme Doppler 
shifts. The performance degradation resulting from using a 
non-uniform window in a situation with little Doppler shift, 
however, is small compared to the performance enhancement 
gained by using a non-uniform window when large Doppler shifts 
exist. In general, in the absence of all information regarding 
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the nature of the received signal, the use of a non-uniform 
window is preferred in all situations. 

For all but the most severe conditions of Doppler shift 
and jamming power, the receiver studied in this thesis 
provides reliable communications through the use of fast 
frequency-hopped spread spectrum. Communications reliability 
does not come without a price. The complexities involved in 
implementing a fast frequency-hopped spread spectrum system 
are enormous and are not worth the benefits in all situations. 
As an example, in order to keep the same number of samples per 
frequency-hop for the DFT, the total symbol duration time 
becomes a linear function of the hop rate. In this thesis the 
hop time is determined to be equal to one second. For the case 
of one frequency-hop per symbol, this translates to a symbol 
duration time of one second, but for the case of L frequency 
hops per symbol, this translates to a symbol duration time of 
L seconds. Thus, a faster frequency-hop rate translates to 
slower, but in many situations, more reliable communications. 

Despite its drawbacks, the use of the DFT based receiver 
may be superior in situations where the content of the message 
is more valuable than the speed of delivery. By employing the 
DFT based receiver with fast-frequency hopping, a system can 
be developed which has the potential to deliver reliable 
communications with considerable anti-jam protection. 
Situations relying on this type of communications protection 
are common in the military environment. 
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In order to gain deeper insight into the relationships 
between the fractional frequency offset, the modulation order, 
the hop rate, and the signal energy-to- jamming noise power 
spectral density ratio, further study should be conducted. 
This study should emphasize the changes in the crossover 
fractional frequency offset observed while varying the key 
parameters listed above. Also, the calculation of the exact 
solution is possible, and should be utilized in any further 
study. Due to limited computational resources, this was not 
possible in this thesis. 

Finally, the work presented here is theoretical in nature 
and should be tested with actual hardware. The assumptions 
made in this analysis were selected to ease the burden of 
computational effort needed to solve the problem. They may not 
correspond to actual conditions needed for optimal 
communications in a real world system. Only building this 
receiver and testing it under realistic conditions will prove 
or disprove the usefulness of the DFT based, noise-normalized 
receiver . 
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